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Abstract. We continue the analysis of the preceding talk with a discussion of the elliptic flow 
and the Hanbury-Brown-Twiss pion correlation radii. It is shown that the thermal model can be 
extended to describe these phenomena. The description of the elliptic flow involves an appropriate 
deformation of the freeze-out hyper-surface and flow velocity. The obtained results agree reasonably 
with the data for soft (< 2 GeV) transverse momenta. For the pionic HBT correlation radii the 
experimental feature that ^ ut/^ s ide — 1 is naturally obtained. The reproduction of individual R side 
and R out can be achieved with the inclusion of the excluded volume corrections, which effectively 
increase the radii by ~ 30%. 



INTRODUCTION 

In the preceding talk [|lj], from now on referred to as (I), it has been shown that the 
thermal approach is successful in the description of the particle ratios and p ± -spectra at 
RHIC. Here we continue our investigation, studying azimuthal asymmetry of the spectra 
and the pionic Hanbury-Brown-Twiss correlation radii. 



ELLIPTIC FLOW 

When the nuclei collide at non-zero impact parameter, b ^ 0, the momentum distribution 
of the produced particles carries azimuthal asymmetry. In general, at mid-rapidity (y = 0) 
we may write the following Fourier decomposition in the azimuthal angle (j>, measured 
from the reaction plane: 



dN 



d 2 p | dy 



dN 



y =o 2np ± dp ± dy 



(l +2v 2 cos20 +2v 4 cos40 + . . .) . (1) 

y=0 



The sines are absent due to the symmetry condition — > — 0, which is simply the 
reflexion with respect to the reaction plane, while the coefficients of cosines with odd 
multiples of vanish for the case of symmetric nuclei and at y = 0, when the symmetry 
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— ► % — holds. The elliptic-flow coefficient, v 2 , can therefore be computed as 



r2x dN 
JO d 2 p ± dy 

V 



2 r2K dN 

JO d 2 p ± dy 



COS 20 d§ 

(2) 



y=0 



The value of v 2 is an important signature of the physics occurring in heavy-ion col- 
lisions. Most importantly, its non-vanishing value indicates that the production mecha- 
nism is not a simple composition of nucleon-nucleon collisions, since in that case the 
asymmetry of production in each such collision would average out practically to zero. 
Thus, interactions with other particles (rescattering, asymmetric collective flow, . . . ) are 
necessary to generate non- vanishing v 2 . In hydrodynamical approaches the elliptic flow 
has been analyzed in many papers, see e.g. [Q, [3J, ||, |5], ^). The coefficient v 2 depends 
on the impact parameter, b, on the transverse momentum p ± , as well as, obviously, on 
the type of the considered particle. All these dependences are measured at RHIC. The 
impact parameter, b, is traded for the experimentally more useful centrality parameter, 
c, which to a very good accuracy is given by [0] 

a2 

(3) 



(2RY 



There are two basic empirical facts from RHIC which we will use in our approach. 
Firstly, v 2 is positive \^ |^, [TI| |TT|, [T2], [T3], [14], [15], P^[F7[], which means that the collective 



flow is faster in the reaction plane than out-of-plane. Secondly, the measurement flT8] ] of 
the azimuthal dependence of the R slde pion HBT radius shows, that the shape of the 
system at freeze-out is elongated out of the reaction plane. We will now use these two 
facts in our choice of the parameterization of the hypersurface at freeze-out. A natural 
extension of Eq. (1.6) is to introduce the azimuthal shape asymmetry, 



r x = Pmax Vl -£ COS(ft, 

Ty = PmaxVl +£sin0, (4) 

with r z and t kept as in the symmetric case of Eq. (1.6) of paper (I). Our convention 
is that r x lies in the reaction plane, and r y is perpendicular to the reaction plane. For 
positive £ this produces elongation out of the reaction plane, as seen in the experiment. 
The four-velocity of Eq. (1.3) is modified as follows: 

u x = Jj r xVl + <5cos0, 

u y = —r y V\ - 8 sin0, 

y N 

1 

u z = — r 7 , 
z N z, 

u t = -t. (5) 
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FIGURE 1. The model fit of the pion, kaon, and proton spectra to the PHENIX data for ^/s NN = 
130 GeV [|l5]| at three centrality bins, arranged top to bottom. Negative (positive) hadron are shown in the 
left (right) side. The optimum values of the size parameters t and p max are given for each centrality bin. 



The normalization N is such that u^u^ = 1. Positive 5 means faster flow in the reaction 
plane, which corresponds to positive v 2 - Certainly, the choice (|],|5|) is by no means 
unique, but it grasps the essential empirical features. 

The modified expansion model has four geometric parameters: T, p mSLX , £ and 8. These 
parameters depend on the centrality parameter, c. Fortunately, the effect of £ and 8 on 
the ^-averaged spectra, dN / (2np ± dp A _dy)\ y _ Q , is negligible and enters at the level of 
a few percent. Thus we may first fit T and p max to the ^-averaged p ± -spectra at various 
centrality parameters, assuming for the moment vanishing e and 8. The result is shown 
in Fig. [I]. We note that the fit works as good as for the most-central case presented in 
(I). The optimum values of parameters are collected in Table [j], where they are also 
compared to the minimum-bias fit and the joint fit to the most central PHENIX [ JT3[ ] and 
STAR [ |T9| ] data. In fact, the qualitative dependence of % and p max on c is as expected: 
the larger c, i.e. the more peripheral collision, the smaller values of the size parameters. 
Figure ^ visualizes this dependence. In Table [T] we also list the ratio of p mSLX / T, and the 
maximum and average values of the flow parameter, /3. Interestingly, these quantities 
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FIGURE 2. The dependence of the size parameters T and p max on the centrality parameter c, as fitted 
to the PHENIX data on ^-integrated /? ± -spectra at v /%^ = 130 GeV Q. 



depend very weakly on c. They are defined as follows: 

1 V^+Pmax' 



/* Pmax r f Pmax 

(jB ± ) = / rdr^==/ rdr. (6) 
io Vr + r Jo 



Ideally, the dependence of the e parameter on c should come from the measurement 
of i? side at various centralities. Then the model evaluation of this quantity would allow to 
fit independently e(c) to the data. We hope to be able to proceed in such a manner in the 
future. Unfortunately, no necessary experimental results are available at the moment. 
In this circumstance we take a reasonable theoretical estimate for e(c) based on Ref. 
[[§, which leads to e = 0.1, 0.21, and 0.35 in the centrality bins - 15%, 15 - 30%, 



TABLE 1. Comparison of the size parameters obtained by fitting 
particle spectra at various values of th centrality parameter c. Their 
ratio, as well as the maximum and average value of the flow parameter, 
/3 ± , are also given. 
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and 30 — 60%, respectively. Finally, the parameter 5 is obtained by fitting the model 
predictions to the v 2 measurements. Our approach includes, as described in detail in 
(I), the decays of resonances. The calculation is straightforward and very similar to the 
one discussed in (I), although takes a much longer computer time due to a lower degree 
of symmetry. The technicalities will be presented elsewhere. The simplified results for 
v 2 presented here include all resonances up to m A = 1.232 GeV, and do not take into 
account three-body decays. 

Figure |3] shows the result of our calculation for three different centrality bins. The 
elliptic flow coefficient grows with the momentum. In continues to grow for large mo- 
menta, where saturation is seen in the experiment, however the thermal model cannot be 
trusted at momenta larger than about, say, 2 GeV, where hard dynamics is important. We 
observe that the effects of resonance decays, large in both the numerator and denomi- 
nator of Eq. (Q), cancel to a large degree in the ratio, and the net effect in v 2 is small. 
The dependence of v 2 on centrality for various transverse-momentum bins is show in 

°- 15 30-60%/ 
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FIGURE 3. Dependence of v 9 (accumulated from all particle species) on the perpendicular momentum, 
p ± , for three centrality bins: 0~- 15% (bottom), 15 - 30% (middle), and 30 - 60% (top). Thick (thin) 
lines correspond to the calculation with (without) resonance decays. Experimental points come from the 
PHENIX collaboration at y/s NN = 130 GeV [[17]]. The taken values for the shape-asymmetry parameter 
£ are, from the lowest to highest centrality bin, 0.1, 0.21, and 0.35, while the fitted values of the flow- 
asymmetry parameter 8 are 0.145, 0.34, and 0.35, respectively. 

Fig. |J We note that our calculation works well except large c and p^, i.e. except for 
high-momentum particles from most peripheral collisions. 

In Fig. H we compare the model predictions for v 2 integrated over c for identified 
particles: pions, kaons, and protons. The general characteristics, with lighter particles 
having larger v 2 , are reproduced. The agreement is not satisfactory only for the case of 
protons at lowest centrality, where the data is compatible with zero. 

To summarize this part we note that 

1 . Elliptic flow can be introduced in the thermal approach by suitably modifying the 
freeze-out hypersurface. 
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FIGURE 4. Dependence of v 2 (accumulated from all particle species) on the centrality parameter for 
various transverse-momentum bins: 0.4 — 0.6 GeV (bottom), 0.6— 1 GeV (middle), and 1 —2.5 GeV (top). 
Crosses show the model calculation, while other symbols are the experimental points from the PHENIX 
collaboration at ^/s NN = 130 GeV [|l7||. The e and 8 parameters are the same as in Fig. ||. 



2. The shape-deformation parameter, e(c), should and hopefully will be taken inde- 
pendently from future data on the azimuthal asymmetry of the R slde HBT correla- 
tion radius at various centralities. 

3. The velocity-asymmetry parameter, 8(c), can be fitted to reproduce v 2 . Predictions 
for identified-particle v 2 follow and are reasonable. 

4. The p ± -dependence shows a monotonic growth, which agrees with the data at lower 
values of p ± , but certainly fails to produce saturation at large momenta, where the 
thermal approach is not applicable. 

5. Resonance decays do not have a very large effect on v 2 . 



HBT RADII 

Now we pass to the description of the pionic Hanbury-Brown-Twiss correlation radii 
(for a review of the problem see, e.g., [|20|]). The studied object is the two-particle 
correlation function for identical particles, in the present case or . It is 

given by 

{n- n~ } 

C{q,P)= r V 2 v (7) 
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FIGURE 5. Dependence of v 2 integrated over c for identified particles on the transverse momentum, 
p ± . Thick (thin) lines correspond to the calculation with (without) resonance decays. Experimental points 
come from the STAR collaboration at ^/s^J =130 GeV [O]. The taken values for the shape-asymmetry 
parameter e are, from the lowest to highest centrality bin, 0.08, 0.25, and 0.52, while the fitted values of 
the flow-asymmetry parameter 8 are 0.15, 0.35, and 0.52, respectively. 




where {.} denotes averaging over events, p x and p 2 are the momenta of the pions, 
q = p 2 —p l , and P = p l + p 2 - We use the Bertch-Pratt parameterization [ |2T| , |Z2| , |Z3| ], 



First, let us briefly recall the experimental highlights. Two facts came as a great 
surprise with the RHIC data. First, the ^ side and R out radii practically do not depend 
on the collision energy [Q, and acquire similar values from AGS to RHIC, despite 
the increase of the energy by almost two orders of magnitude. Secondly, the ratio of R out 
to R slde , which can be interpreted as a measure of the duration time of the freeze-out, 
is close or even less than one. This is in contradiction to anticipations from numerous 
hydrodynamic simulations, which had predicted ^ ut/^ s ide significantly larger than one. 



Our model evaluation of the correlation function is performed according to the for- 
malism of Ref. with a technical approximation of neglecting the finite life-time of 
resonances. In that case 



_ \JdL(x)-u(x)e i i- x S(P-u(x))\ 2 
C{q,P) = 1 + j dL . uS ((p + zy u(x))JdZ-uS((P-^) ■ u(x)) ' (9) 

where the source function is 

S(p-u) = t— -t T r' p ' M_ ' 1 ^ 7 ' + contribution from resonances. (10) 

As discussed in Ref. [[251], our model values for the geometric parameters % and p max 
are low, of the order of the size of the colliding nuclei. As a result, the values of the 
i? side and i? out HBT radii obtained with the procedure described above are about 30% 
too low compared to the experiment. The problem can be alleviated with the inclusion 
of the excluded- volume (Van der Waals) corrections. Such effects have been realized 
to be important since the early thermal studies of the particle production in relativistic 
heavy-ion collisions [ [261 , WA 0], where they led to a significant dilution of system. 
In the case of the Boltzmann statistics, which is a very good approximation in our case 



[pOQ, the excluded volume corrections bring in a factor p8fl 



(11) 



I+Ljvje-Wn/ 

into the phase-space integrals, where P is the pressure, v- = 4|^r? 
volume for the particle of species i, and n ; is the density of particles of species i. The 
pressure can be calculated self-consistently from the equation 

P = %P?{T^ i -Pv i /T) =^Pf{T^)e- pv >l\ (12) 

i i 

where Pf denotes the partial pressure of the ideal gas of hadrons of species i. With the 
simplest assumption that the excluded volumes for all particles are equal, r- = r, the 
excluded-volume correction manifests itself as a common scale factor, which we may 
denote by S~ 3 . The Frye-Cooper formula can then be written in the form pi], [32|1 



dN, , [+°° /-Pmax/T 

—5 — — = t / da,, / sinha , d sinha ,) 

d z p ± dy J-™ ll Jo 

x / dZp-uS-iffip-u), (13) 
Jo 

where p ■ u = m ± cosha||Cosha ± — p ± cos^ sinha^. As can be immediately seen from 

this expression, the presence of the factor S~ 3 in Eq. Q) may be compensated by 
rescaling p and % by the factor S. That way the system becomes more dilute and larger 
in such a way, that the particle multiplicities and the spectra are left intact. 
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FIGURE 6. Model predictions for the pionic R side and R out HBT correlation radii (top two panels), 
and their ratio (bottom panel), confronted with the PHENIX data Au+Au data at ^/s^ =130 GeV and 
average centrality 10%. The quantity k ± is the total momentum of the pion pair. 




With our values of the thermodynamic parameters (T) jU-) = 80 MeV/fm 3 , and 
we find S = 1.3 with r = 0.6 fm. Such a value of the excluded volume is compatible 
with values typically obtained in other calculations. The increase of the size parameters 
by 30% is what we need to bring the values of R slde and R out up to the experimental ball 
park. 

Our results are shown in Fig. We note that the agreement with data is very rea- 
sonable. In particular, the ratio of ^ ut/^ s ide 1S c l° se t0 one > an d drops below one at 
larger values of the pair momentum k ± . The plots of R sids and R out include the excluded- 
volume correction factor S = 1.3. We observe that these radii decrease with the pair 
momentum k ± , although somewhat slower than indicated by the data. We note that the 
radius ^i ong cannot be reliably evaluated in our model. This is due to the assumption of 
the boost invariance, cf. Eq. (1.2), which leads to too large values of /? long . 



CONCLUSION 



To conclude, we list the main results of our approach: 

1. The thermal model works for the particle ratios, see (I). 

2. Supplied with expansion, it works for the p ± - spectra, see (I). The complete treat- 
ment of resonances is essential, and the assumption of single freeze-out leads to 
good predictions. Moreover, the strange particles, including £1, are described prop- 
erly with no need for extra parameters. 

3. Supplied with azimuthal asymmetry, the model can be used to describe the elliptic 
flow at moderate transverse momenta, up to p ± ~ 2 GeV. The fitted values of the 
velocity asymmetry parameter, 5, are reasonable. 

4. Supplied with the excluded- volume corrections, the model works also for the HBT 
radii i? side and i? out . In particular, the ratio of ^ ut/^ s ide * s c l° se to 1. 

5. The description is efficient, involving two thermal parameters T and jj. g , two size 
parameters % and p max , and in the case of azimuthal asymmetry, two deformation 
parameters, e and 8. 

6. Finally, we note that the model also works for the case of RHIC at *Js^ = 



200 GeV A as well as for SPS at v ^ = 17 GeV Q3J, 
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